Topological Quantum Computation by Manipulating Quantum Tunneling Effect of the 

Toric Codes 
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Quantum computers are predicted to utilize quantum states to perform memory and to process 
tasks far faster than those of conventional classical computers. In this paper we show a new road 
towards building fault tolerance quantum computer by tuning quantum tunneling effect of the 
degenerate quantum states in topological order, instead of by braiding anyons. Using a designer 
Hamiltonian - the Wen-Plaquette model as an example, we study its quantum tunneling effect of 
the toric codes and show how to control the toric code to realize topological quantum computation 
(TQC). In particular, we give a proposal to the measurement of TQC. In the end the realization of 
the Wen-Plaquette model in cold atoms is discussed. 
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Quantum computers are predicted to utilize quantum 
states to perform memory and to process tasks far faster 
than those of conventional classical computers. Various 
designs have been proposed to establish a quantum com- 
puter, including manipulating electrons in a quantum dot 
or phonon by using ion traps, cavity QED, or nuclear spin 
by NMR techniques. Recently, people find that it may 
be possible to incorporate intrinsic fault tolerance into a 
quantum computer - topological quantum computation 
(TQC) which has the debilitating effects of decoherence 
and free from errors. The key point is to store and ma- 
nipulate quantum information in a "non-local" way, as 
means that the "non-local" properties of a quantum sys- 
tem remain unchanged when one does local operations on 
it. An interesting idea to realize fault-tolerant quantum 
computation is anyon-braiding proposed by Kitaev [l|, Q. 
He pointed out that the degenerate ground states of a 
topological order make up a protected code subspace (the 
topological qubit) free from errors^ 0]. 

It is known that there are two types of topological or- 
ders in two dimensions spin models - non-Abelian topo- 
logically ordered state and Z2 topologically ordered state. 
Z2 topological order is the simplest topologically ordered 
state with three types of quasiparticles: Z2 charge, Z2 
vortex, and fcrmions[15|. Z2 charge and Z2 vortex are 
all bosons with mutual ir statistics between them. The 
fcrmions can be regarded as bound states of a Z2 charge 
and a Z 2 vortex. Recently, several two dimensional ex- 
actly solved spin models with Z 2 topological orders were 
found, such as the toric-code model the Wen's pla- 
quette model 
eycomb lattice 




1J] and the Kitaev model on a hon- 
One can operate on the protected 
code subspace of Z2 topological order by creating an ex- 
citation pair, moving one of the excitations around the 
torus, and annihilating it with the other. However, the 
non-local operations do not form a complete basis. 

On the other hand, in non-Abelian topological orders 
the elementary excitation becomes non-Abelian anyon 
with nontrivial statistics. Now people focus on real- 
izing TQC by braiding non-Abelian anyons. The de- 



generate states undergo a nontrivial unitary transforma- 
tion when a non-Abelian anyon moves around the other. 
One can initialize, manipulate and measure the degener- 
ate ground states with several non-Abelian anyons]^, @1 ■ 
Along the road, the key point is to manipulate sin- 
gle quasi-particle which becomes a hot issue recently 

In particular, 
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I12L 1 1 31 ] . In particular, the idea of 
creation and manipulation of anyons based on controlled 
string operators is proposed [il [Iol|. 

Kitaev have noted the degenerate ground states (on a 
torus) as the toric code. In this paper We will design 
TQC by manipulating the degenerate ground states of 
Z2 topological orders (the toric code), instead of that by 
braiding anyons in a non-Abelian topological order. The 
key point to manipulate the degenerate ground states is 
to tune quantum tunneling effect by controlling external 
field on spin models. Firstly the effective theory of the 
toric code in Z2 topological orders is formalized. Sec- 
ondly, by using the Wen-Plaquette model as an example, 
we show how to control the toric code by tuning the tun- 
neling of the degenerate ground states [3, [lj, [l_5j]. In this 
part, we will concentrate on the measurement of toric 
code. Finally we give a short discussion on the realiza- 
tion of the Wen-Plaquette model in cold atoms. 

The effective model of the toric code : For Z2 topolog- 
ical orders on an even-by-even (e * e) lattice on a torus, 
the ground state degeneracy is 40, [l4, 15 1. The four de- 
generate ground states are denoted by | m, n), m = 0, 1 
and n = 0, 1. m, n) have different boundary conditions 
of fermion's wave- functions ip(x, y) as 

1>{x, y) = (-l) m *P(x, y+L y ), $(x, y) = (-l) n ^{x+L x , y). 

(1) 

Physically, topological degeneracy arises from presence or 
the absence of ir flux of fermions through the two holes 
of the torus. The degenerate ground states make up two 
qubits which can be mapped onto quantum states of two 
pseudo-spins fi and f 2 , | 0, 0) — >-|t)i® I T> 2 , | 1, 0) — >\i 

>i®|T>a, |o, 1) -|T>i® U> 2 , 1 1, i> ->|i>i® |i> 2 . 

It is known that the degenerate ground states with 
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topological orders have same energy in thermodynamic 
limit. In a finite system, the degeneracy of the ground 
states is (partially) removed due to tunneling processes, 
of which a virtual quasi-particle moves around the torus 
before annihilated with the other onep], 0, Q. In Z% 
topological orders, there are nine tunneling processes de- 
noted by Q+y, C x fl Cj, C x f +y , Cjf, C x c +y. : that 
correspond to virtual ^-vortex, fermion, Z2 charge prop- 
agating along e x , e y and e x + e v (or e x — e y ) directions 
around the torus, respectively. For example, the process 
C x becomes the unitary operation as 

(lt)i, IDi, lt> 2 , U> 2 )-»(U>i, IT)i, lt> 2 , |l) 2 ). 

(2) 

1 



Here J XXl J yy , J zz , J zx , J xz , hf, h x . hf, h\ are deter- 
mined by the energy splitting of the degenerate ground 
states from the nine tunneling processes in a Z2 topolog- 
ical order. 

Manipulating the toric code by controlling tunneling 
splitting : To design a topological quantum computer, 
one needs to do arbitrary unitary operations on the toric 
code. To emphasize this point, we introduce a con- 
cept 'Controllable Topological Order (CTO)\ In a CTO, 
quasi-particles' dispersions and the energy splitting of the 
degenerate ground states can be manipulated [9I llO] . 

In the following part, based on an example of control- 
lable Z2 topological order - the Wen-plaquette model, 
we demonstrate how to do TQC. The Hamiltonian of the 
Wen-plaquette model is0, [3, EH 

H Q = -gY i F i . (4) 

i 

Here F t = a x a y + ^a x + ^ +iy a y + ~ y and g > 0. a, are Pauli 
matrices on sites, i. The ground state denoted by F{ = +1 
at each site is a Z2 topological state with the ground 
state e nerg y E g = —gN where N is the total lattice 
number [la]. On an e * e lattice, the four degenerate 
ground states make up two qubits. However, on even- 
by-odd (e * o) , odd-by-even (o * e) and odd- by-odd (o * o) 
lattices, the ground state degeneracy is two instead of 
four. 

In this model Z2 charge is defined as Fi^ even = — 1 at 
even sub-plaquette and Z2 vortex is Fj eo dd = — 1 at odd 
sub-plaquette. The fermions can be regarded as bound 
states of a Z2 charge and a Z2 vortex on two neighbor 
plaquettes. These quasi-particles in such exactly solved 
model have flat band. In other words, the quasiparticles 



Hence we can use the a pseudo-spin operator rf ® 1 to 
represent C x . Similarly Q, C x+ v , C x , C%, C x+ y , C x f , Cj, 
C x+y can be denoted by 1 ® rf , rf <g rf , rf ® rf , rf g) rf , 
rf ®rf , l®rf , rf £3>1, rf (8rf , respectively. Among them, 
one can choose four basic processes, C x , C%, C5, C y , and 
check the ground state degeneracy from the commutation 
relations between C x and C y , that obeys the Heisenberg 
algebra, C x C y = e l7r C y f C x . 



Hence the effective Hamiltonian of the toric code is 
obtained in term of the pseudo-spin operators : 



(3) 

I 

cannot move at all. Under the perturbation 

H' = h'Y^af + fc*X> i *, (5) 

i i 

the quasiparticles begin to hop. The term h x ^2 af drives 

i 

the Z2 vortex, Z2 charge and fermion hopping along di- 
agonal direction e x + e y . The term h z ^ o~f drive fermion 

i 

hopping along e x and e y directions without affecting Z2 
vortex and Z2 charge. One can see the detailed descrip- 
tion of the three kinds of quasi-particles in Ref . [l5| • 

As a result, there exist five tunneling processes under 
the perturbation H', C x+y , C x f +V , C], C y , C x+y . Let us 
calculate the ground state energy splitting from a higher 
order (degenerate) perturbation approach. From the tun- 
neling processes of Z 2 vortex and Z 2 charge [USQ, C x+y 
and C x+y , one can determine J xx = J yy — J ~ g(^y) L ■ 

For L y lattice on a torus, L is equal to " f y 

where £ is the maximum common divisor for L x and 
L y . From the tunneling process of fermion, C x , C y and 

C x+V , one can obtain hf ~ g(^) Ly , K ~ 5(ff) Lx and 
J zz ~ g(^) 2L + 9 (|j) L x+£ B _ Other parameters are all 

zero, J zx = J xz = hf = /if = 0. Then the effective model 
is simplified into 

H cfi ~ J(rf-rf + rf-rf) + J zz rl-rl + +^rf + ^r|. (6) 

On e * o, o * e or 0*0 lattices, the effective model 7Y e ff 
becomes more simple. For example, the two degenerate 
ground states on an L*L ( L is an odd number) lattice are 
|t)i® |T)2 and ID2 which are denoted as ||) and |J.) 

in the following parts of this paper. And the pseudo-spin 
operator from different tunneling processes is denoted by 
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FIG. 1: The energy splitting between the two degenerate 
ground states SE from the exact diagonalization numerical 
results (g is set to be unit, g=l). Here N * M denotes a 
N x M lattice. 

t. Then the dynamics of a single qubit can be described 
by a simple effective pseudo-spin Hamiltonian 

Weff = J(r x +Ty) + J zz T z . (7) 

Fig.l shows the energy splitting (scaled by g ) of the two 
degenerate ground states from the exact diagonalization 
numerical results. 

Then by adding the specific perturbations to the Wcn- 
plaquettc model, H' , one can change different quasi- 
particles' hopping and then manipulate the toric code 
by controlling tunneling splitting of degenerate ground 
states. Such dramatic property in the Wen-plaquette 
model gives an example of so-called controllable topolog- 
ical orders. Similar properties have been used to creat 
and manipulate anyons in the Kitaev modeljij. Iioj. 

Topological quantum computation : In this section we 
focus on two degenerate ground states on an o * o lat- 
tice (|t) and |J,)) and show the initialization, the unitary 
transformation and the measurement. 

Firstly we initialize the system into the quantum state 
||). This process will occur according to the Hamilto- 
nian H'= h{t)Y,Vz where h(t) = e~ t/to - 1. At first 

i 

t — ► — oo, the ground state is the spin-polarized state 
of a z eigenstates \ipo). Since the effective Hamiltonian of 
the toric code is TC c s = h z ■ t z , the time-evolution oper- 
ator U{t) = e s in the topologically ordered phase 
becomes a projection operator of the pseudo-spin as 
U(t)\ip) —►|T)- The system will eventually evolve adia- 
batically and continuously from \ipo) into the final state 
U(t = O)|-0o) —►!!)• It) becomes the initial state prepared 
for TQC. 

Secondly we show unitary operations of the TQC [1,0, 
[Tlj . A general pseudo-spin rotation operator is defined 
as 

U e , v = er^e-rM-'+^e-^ ( 8 ) 



where 7 = J zz At 7 , 6 = J zz Atg and ip = JAt v . For 
example, the Hadamard gate can be a special pseudo-spin 
rotation operator, Ug tV {"f= j, 9= ip = j). To design 
the Hadamard gate, we may apply the external field along 
z-direction at an interval Atg = -^j— firstly. Then, we 
swerve the external field along x-direction at an interval 
At v = jj. Finally, the external field along z-direction is 
added at an interval Ai 7 = jj—- Using similar method, 
one can reach certain quantum operations demanded by 
TQC and have the ability to carry out arbitrary gate 
onto the toric code at will, a ||) + (3e % ^ ||) with a, (3 > 
(a 2 +/3 2 = l). 

Thirdly let us discuss the measurement of an arbi- 
trary state I vac) = a |T) + Pe^ \[). The interference 
from Aharonov-Bohm (AB) effect allows one to design an 
experimentally observable distinction between the pro- 
cesses with or without a 7r-fiux inside the loop. To deter- 
mine a, (3 and </>, we need to observe both fermion inter- 
ference and Z2 vortex interference. Fig. 2 is a scheme to 
show the AB interference on a torus. 

Let's explain how to determine a and (3 by AB effect 
from fermion-interference. To observe the AB interfer- 
ence, we add a small external field h x — ► and h z = 0. 
Now Z2 vortex, Z2 charge and fermion hop along diagonal 
direction e x + e y . On a torus, there exist two symmetrical 
paths between two sites ( i and j ) on opposite positions 
on a torus, 71 and 72. Then the two trajectories will 
contribute to the transition amplitude Tij according to : 

r *-j = l*«r + l*35l 2 + 2e l*S*«l ( 9 ) 

where and ^J^- are the wave functions of fermions of 
the two trajectories. For the ground state |T), e is unit, 
e = 1. However, for the ground state ||), e = —1. Then 
we can distinguish these two cases. For two symmetrical 
paths = ^7j = */' we § et a probability a 2 for |T) 
with Tij = At 2 and a probability (3 2 for \l) with T itj = 0. 

On the other hand, one can determine the parameter 
4> by observing Z 2 vortex interference from i to site j. 
The wave function of Z2 vortex has a periodic boundary 
condition along x direction for the ground state |T') = 
-j= |T)+ |l) and an anti-periodic boundary condition for 
the ground state |J/) = |T)— ||). Then an arbitrary 
state I vac) = a |T) + (3e 1 ^ ||) is re-written into 

I vac) = ^l+apcosfa**' IT') + ^ - a/? cos ^e^" ||') 

(10) 

where </>> = arctan( ? ^_) and 0" = arctan(^^). 
For two symmetrical paths = ^Tj = ^v, we get a 
probability (| + a/3 cos (f)) for |T') with Tij = At 2 and a 
probability ^ — a/3 cos </) for |T') with T; j =0. As a result, 
we determine the parameters a, (3 and <j) of an arbitrary 
state I vac) = a \{) + f3e^ ||). 

For Z2 topological orders on an e * e lattice, the pro- 
tected subspacc becomes two qubits. To do arbitrary 
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FIG. 2: The scheme of the interference on a torus. 

unitary transformation on the protected subspace, one 
needs to apply the external field on a sub-lattice, 



H' 



i£odd 



+ hJ2 ot + h'^a?. (11) 



zGodd 



This is because the perturbation h °~i + h X) °\ 

iEodd iGodd 

drives only Z2 vortex hopping without affecting Z% 
charge. Hence the nine parameters of the effective 



h\L 
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J ZX ~ 

h\L v 



pseudo-spin model become J xx 

J XZ = 0, J zz = g(^) L * +L \ hi - yy-j „ 2 

hf = g(j-) Ly and h\ = g(^-) Lx ■ Furthermore, on a man- 
ifold with higher genus, the degenerate ground states 
become multi-qubit, which can be mapped onto multi- 
pseudo-spin model. By the same method one can do 
TQC on the multi-qubit. 

Realization of the Wen-plaquette model : Finally we 
discuss the realization of the Wcn-plaqucttc model in a 
optical lattice of cold atoms. Because the Wen-plaquette 
model can be regarded as an effective model of the Kitaev 
model on a two dimensional hexagonal lattice, one may 
realize the Kitaev model firstly. The Hamiltonian of the 
Kitaev model is 

h = Yl ( J *°i,i°i+i,i + vj-i^j.i + J«°i,i°i,j+i) 

even 

(12) 

where j and I denote the column and row indices of the 
lattice. In the limit J^^- 3 Z ~ J y of Kitaev model, the 
effective Hamiltonian of the Kitaev model is simplified 
into the Wen-plaquette model 



Ha 



t2t2 

— 3 Z^ CT lcft(i 



16 J 



x y y 

<T right(i) Cr up(i) C7 down(i) ■ 



(13) 



Then one can use the Kitaev model in the limit 3 Z ~ 
J y on a torus to do the TQC. The realization of the 
Kitaev model on a two dimensional hexagonal lattice 
has been proposed in Ref.0, [l(| EH 0, 17 [■ The es- 
sential idea realizing the Kitaev model is to induce and 
control virtual spin-dependent tunneling between neigh- 
boring atoms in the lattice that results in a controllable 
Heisenbcrg exchange interaction. 



Summary and discussion: By controlling the tunneling 
processes, one can do unitary transformation on a topo- 
logical protected qubit (the toric code), as paves a new 
road towards TQC rather than by anyon-braiding. By 
using a designer Hamiltonian - the Wen-Plaquette model 
as an example, we show how to control the toric code 
to realize TQC. In particular, we give a proposal to the 
measurement. 

In the end we give a short comment on the advan- 
tage and the disadvantage of TQC by tuning quantum 
tunneling effect. For TQC by anyon-braiding, one needs 
to manipulate single quasi-particle which demands new 
techniques. On the contrary, for the TQC by tuning tun- 
neling, one needs only adjust the global field strength (or 
direction). However, it is a true challenge to realize the 
designed spin model on a manifold of higher genus in the 
optical lattice of cold atoms. Such unsolved issue will be 
worked out in the future. 
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